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We realize a magnetoresistive emulsion analyzer capable of detection, multiparametric analysis and sorting 
of ferrofiuid- containing nanoliter-droplets. The operation of the device in a cytometric mode provides high 
throughput and quantitative information about the dimensions and magnetic content of the emulsion. Our 
method offers important complementarity to conventional optical approaches involving ferrofluids, and 
paves the way to the development of novel compact tools for diagnostics and nanomedicine including drug 
design and screening. 

A substantial part of nanomedicine focuses on drug design and employing magnetic nanomaterials in 
combination with multifunctional polymers, lipids or proteins^'^. For example, ferrofluid-based liposomes 
are capable to respond to external stimuli like magnetic fields, which allows the remote control over the 
chemical species^ Ferrofluids are widely used for chemotherapies^, in vivo magnetic resonance imaging (MRI)^ 
and as labels in immunological tests, i.e. magnetic enzyme-linked immunosorbent assays (ELISAs)'^. 

High-throughput drug discovery requires a tool which is able to analyze, manipulate, and sort objects, e.g. 
liposome capsules or emulsions containing known doses of magnetic nanoparticles associated with medications. 
A promising route relies on the implementation of magnetically-labeled biochemical species in combination with 
magnetoresistive sensors. In-flow detection based on giant magnetoresistive (GMR)^ or tunneling magnetore- 
sistive (TMR) sensors^ has been demonstrated recently, representing first steps towards the integration of 
magnetic sensor elements into fluidic channels. 

However, up to date, progress is restrained to the mere sensing and counting of magnetic objects. Advanced 
and quantitative analysis, which is the main advantage of e.g. optical flow cytometry^, has not been explored so far 
with magnetic sensors. Apart from detection and analysis, sorting of species is invaluable in diagnostic devices. 
The sorting of biochemical objects has been dominated by conventional fluorescence-activated cell sorting 
(FAGS) techniques^'. 

Here and for the first time, a magnetic emulsion analyzer capable of multiparametric studies and sorting is 
realized. Ferrofiuid emulsion droplets of different size and concentrations of magnetic nanoparticles are pro- 
duced. The droplet volumes range from 20 nl up to 400 nl - ideally suited for culturing monoclonal cell 
populations^^ and monitoring bacteria growth^^. Multiparametric analysis as well as magnetic sorting of ferrofiuid 
droplets of different size are performed to demonstrate the potential of OMR sensors for magneto -cytometry 
studies. 



Results 

The concept of the magnetoresistive emulsion analyzer capable of multiparametric studies and sorting of ferro- 
fiuid droplets is presented in Figure la. Water-in-oil emulsion droplets of various sizes with encapsulated 
ferrofluids are produced in a T-junction geometry. An example of the flow profiles of the two immiscible phases 
is displayed in Supplementary Figure 6. 

The key component in the experimental setup is the GMR-based sensing platform. The sensor is composed of a 
Py/Cu (Py = NigiFcig) layer stack coupled at the 2nd antiferromagnetic maximum (AFM). The magnetoresistive 
characteristics (OMR curve and sensitivity) are shown in Figure lb. A GMR ratio of 1 1% and sensitivity of up to 
26 T"Mn a magnetic field of 0.9 mT is achieved. The high sensitivity of the GMR sensor and the hysteresis -free 
sensor response under cycling magnetic field are crucial for the magneto -cytometric detection. The detection is 
based on a Wheatstone bridge geometry (Supplementary Figure 2). A lock- in amplifier is used to improve the 
signal-to-noise ratio. Affiliated to the sensing platform is a digital camera to monitor the flowing droplets. 
Recorded videos and images are correlated with the voltage signals obtained with the GMR sensor. The 
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Figure 1 | (a) Concept of GMR-based droplet analyzer for multiparametric analysis and sorting, (b) Magnetoresistive characterization of GMR sensor. 
Inset: Photograph of the patterned sensor prepared for measurements in 4-point configuration, (c) Real time detection of a train containing small 
and large ferrofluid emulsion droplets. Sorting of small (dl) and large (d2) droplets at the T-junction. 



water-based ferrofluid (EMG 700, Ferrotec) does not retain reman- 
ence without applying an external magnetic field. Therefore, a bias- 
ing field is applied to achieve a net magnetization of the ferrofluid, 
which also serves to operate the sensor in its most sensitive range. 

First, we demonstrate the capability of the GMR-based emulsion 
analyzer to detect and examine ferrofluid droplets of different sizes 
and concentrations. Each detection event is represented by a peak 
(Fig. Ic; Supplementary Figure la) characterized by its amplitude 
and the full- width at half maximum (FWHM). These characteristics 
serve as fingerprints for a ferrofluid droplet, flow conditions and 
emulsification procedure. For instance, in the case of a fixed sen- 
sor-droplet separation, the amplitude is related to the size and con- 
centration of the magnetic material in a droplet, whereas the FWHM 
contains information about the spread in speed and size of a droplet. 
Since the following experiments were performed at a fixed flow rate, 
modifications of the FWHM are mainly attributed to the variation in 
droplet size. 

Discussion 

The multiparametric density plots produced for different volumes of 
ferrofiuid droplets (concentration: 15 mg/ml) are summarized in 
Figure 2a. The broad distribution of the FWHM and amplitude for 
the smallest droplets (black cloud. Fig. 2a) are consistent with the real 
time detection results (Supplementary Figure la). The correspond- 
ing broad standard deviations of the fitted histograms in Figure 2c 
reflect predominantly the clear polydispersity of droplet sizes, and is 
to lesser extent ascribed to the inhomogeneity of the ferrofluid con- 
centration between droplets. The polydispersity is partly related to 
the spontaneous fusion of droplets smaller than the channel dimen- 
sion. With increasing size (compare red and green clouds), the drop- 
lets completely flll the channel cross-section resulting in a more 
uniform droplet train as reflected by the narrower distribution of 
the scattered data points around the cloud center (Supplementary 
Figure 2i). Once the droplet is longer than the width of the magnetic 
sensor, the sensor voltage change (amplitude) reaches its saturation 
value and becomes insensitive to the size of the droplets. However, 



the information about the size of droplets can still be extracted from 
the FWHM (compare green and blue clouds. Figure 2a). An increase 
of the FWHM magnitude is observed when the droplet volume 
increases from 120 nl to 270 nl. 

Multiparametric analysis performed for a train of droplets 
(volume: 160 nl) with different concentrations of ferrofluid is shown 
in Figure 2b. The clouds of data points for different concentrations 
naturally fall into the same FWHM region, and the standard devi- 
ation of the amplitude and FWHM levels off (Supplementary Figure 
2j), which implies a relatively high homogeneity of the droplet size 
and concentration of the ferrofluid for each experiment. The FWHM 
magnitude within the cloud is mainly determined by the deviation of 
droplets speed. In contrast, due to the scaling of the voltage change 
with the magnetic moment, the center of the data point cloud is 
clearly shifted to higher voltage amplitudes for droplets with larger 
concentration of magnetic nanoparticles (Figure 2d). 

The magnetic droplet analyzer can be extended towards sorting 
emulsions with respect to a specific property addressed in a multi- 
parametric analysis. The sorting module consists of a split junc- 
tion and a solenoid 3 -way isolation valve installed behind the 
GMR sensor. For demonstration purposes, we performed sorting 
of droplets of different sizes, which provide various readouts of 
the GMR sensor (Figure Ic, Supplementary Movie 1). Small and 
big droplets were produced in a continuous mineral oil flow (Qon 
= 200 nl/s) using two separated T-junctions (Supplementary 
Figure 5b) Ferrofluids with a concentration of 75 mg/ml were 
used for the sorting experiments. The small droplets (volume: 
20 nl) cause an average voltage change of about 7 jiV, while the 
big droplets (volume: 500 nl) provide an average positive voltage 
change of about 30 |iV. Droplets of different size can be sorted 
into two separate channels via controlling the isolation valve to be 
open only when big droplets were passing by the GMR sensor 
(Supplementary Movie 2). It is foreseeable that various types of 
emulsion droplets can be sorted with this tool by implementing 
additional sorting criteria, e.g. dose of the encapsulated magnetic 
nanoparticles. The demonstration of magnetic -activated sorting 
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Figure 2 | Multiparametric density plots for the detection of magnetic droplets which are produced with (a) different volume of ferrofluid droplets 
(concentration: 15 mg/ml) and (b) different concentrations of ferrofluid in droplets (volume: 160 nl). The flow rate of mineral oil is Qoh = 500 nl/s. 
Circles are guide to the eyes. Histograms of the voltage amplitude change with (c) different volume of ferrofluid droplets and (d) different concentrations 
of ferrofluid in droplets. These plots collect measurement data of about 100 droplets. 



opens up promising perspectives for drug design and screening, 
which could be beneficial to the pharmaceutical industry. 

Methods 

Droplet generation module. The ferrofluid droplets are formed on a T-junction 
connector, where aqueous solution of ferrofluid (water-based ferrofluid, Ferrotec 
EMG Series 700) is dispersed in the mineral oil phase (mineral oil, Sigma- Aldrich 
M8410). The droplet generator module (Supplementary Fig. 5a) was assembled from 
several components: the pulse free syringe pumps (neMESYS, Cetoni) for fluid 
injection; the T-junction for droplet formation; interconnections, i.e. flangeless 
fittings (standard %-28 thread), syringes (1000 Series, Hamilton, syringe volume: 
1.0 ml and 2.5 ml); flow circuits that were realized with polytetrafluoroethylene 
(PTFE) tubing with outer diameter (OD) of 0.9 mm and inner diameter (ID) of 
0.4 mm. We used the pulse flow approach to generate droplets of different size for 
multiparametric analysis. An example of the flow profiles of two immiscible phases 
are displayed in Supplementary Figure 6. The period of flow pulses is 7 s, while each 
flow pulse lasts for 500 ms. 

Magnetic detection module. The layer stack of the GMR multilayers was: [Py 
(1.5 nm)/Cu (2.2 nm)]3o/Py (1.5 nm), which was prepared by magnetron sputtering 
technique under high vacuum condition (base pressure: 2.5 X 10"^ mbar). Ar was 
employed as the sputter gas, the pressure of which was controlled at 1.6 X 10"^ mbar. 
Before depositing the layer stack, a thermally oxidized silicon wafer with 1000 nm 
SiOx layer was patterned into rectangular stripe (width: 1 mm, length: 16 mm) by 
standard optical lithography using image reversal technique. The optical lithography 
was performed in the clean room with standard class 100 conditions. Firstly, the 
silicon substrate covered with 1 |j,m oxide was spin coated with Tl-Prime 
(MicroChem Corp.) at 4500 rpm for 35 s, followed by a baking at 120°C for 3 min. 
After that, negative photoresist (AZ5214E) was spin coated onto the substrate at 
4500 rpm for 35 s, which was subsequently baked at 90°C for 4 min. Then the sample 
was exposed with the stripe pattern using MA56 mask (Karl Siiss) aligner for 2 s, and 
post-baked at 120° C for 2 min. Subsequently, the sample was exposed without mask 
for 30 s to realize image reversal. At the end, the sample was developed by AZ726MIF 
developer (Microchem) for 1 min, and rinsed by DI water. 

Magneto -electrical characterization was performed based on conventional 4-point 
method. The GMR sensor was placed in between pole shoes of the electromagnets, 
where uniform in plane magnetic field was cycled between ±25 mT and the sensor 
resistance was recorded. The GMR sensor is a proximity sensor, which responses to 
local magnetic stray fields from the objects. Therefore, a strategy is taken to enhance 
the magnetic stray fields to be detected. As the ferrofluid (EMG 700, Ferrotec) 
employed for present study is a suspension of magnetic colloidal nanoparticles in 
water, it does not retain remanence without applying external magnetic field. Thus, 
an external field should be applied to achieve a net magnetization of the ferrofluid. In 
addition, the external magnetic field can be simultaneously used to bias the sensor to 



the sensitive region. Based on the above considerations, a permanent magnet (AlNiCo 
500, type A 1560, IBSMagnet, length: 60 mm, diameter: 15 mm) was placed below the 
GMR sensor which was mounted on a movable stage. With adjustment of the position 
of the permanent magnet, the sensor was biased to the most sensitive region and the 
stray fields from the magnetized ferrofluid were sufficient to be detected by the GMR 
sensor. Further, the tubing was placed in close contact to the sensor surface to 
enhance the stray fields that can be detected. As the GMR sensor is the proximity 
sensor, the transition from the dynamic to saturation regime (Supplementary Fig. 
lb) is related to the dimension of the GMR sensor but not entirely to the geometrical 
parameters of the tube. Saturation regime is observed when the sensor width 
(dimension of the sensor in the direction of the droplet travel) is equal to the length of 
the droplet. In our experiments, width of the sensor is 1 mm. Therefore, transition to 
saturation is observed when the droplet length is also about 1 mm (Supplementary 
Fig. lb). In this case, the critical volume of the droplet assuming its shape of cylinder 
with two half- spheres at both ends can be estimated by taking into account width of 
the sensor (W) and radius of the fluidic tube (R): Vc = 4/3 7i R' + Pi R^ (W - 2R). For 
the geometry of our experiment, Vc ~ 110 nl. 

The GMR sensor was soldered and connected with 3 other trimmers (Bourns, 
3006P-l-102LFto realize a Wheatstone bridge geometry (Supplementary Fig. 5c). By 
balancing the bridge, the background voltage from the bridge can be minimized to 
enhance the measuring sensitivity. The lock-in amplifier (SRS830) was used to 
amplify the signal and reduce the noise. The whole bridge circuit was powered also by 
lock-in with constant ac voltage of 0.800 V. The modulation frequency was 1 kHz 
and the time constant was 10 ms. The sensitivity of lock- in was set at 500 |iV with a 
resolution of 10 nV, which was sufficient to resolve the voltage signal. The analog 
output from the lock- in was sampled at 512 Hz and transferred to the computer. 
Affiliated to the sensing platform a digital camera (DNT Microdigit) was installed 
which was used to monitor the flowing droplets. Videos and pictures of flowing 
droplets can be correlated with the voltage signals obtained with the GMR sensor. 

GMR-based sorting module. The droplet generation module was refined to form a 
train of droplets of multiple sizes for the sorting (Supplementary Fig. 5a). The 
module was composed of two T-junctions (Gland C2) with different thru-hole 
diameters (CI: 0.5 mm, C2: 0.25 mm). Due to the differences in the size of thru-hole 
diameter, droplets of distinct size can be produced at these two junctions. Big droplets 
were formed in the junction CI using syringe A, while smaller droplets were injected 
into the same fluidic channel through the junction C2 through syringe C. Such fluidic 
assembly resulted in the sequence of the droplets with alternating size, which was 
further used for sorting experiments. The sorting part was composed of a T-junction 
(Supplementary Fig. 5d, thru-hole diameter: 0.25 mm) flanked by two separate 
fluidic channels that were connected to a solenoid 3 -way isolation valve (Amico 
Scientific) as mentioned before. The 3-way valve was with one-way open and one-way 
closed without supplying power. The switching of the valve was controlled by an 
external power source (Agilent Precision Source/Measurement Units B2902A) with 
working voltage of 24 V. In the sorting experiment, applied voltage of 20 V with 
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compliance of 0.5 mA was sufficient to switch on the valve. The algorithm of 
magnetic sorting is described as follows: ferrofluid droplets of different size are 
detected by the GMR-based millifluidic droplet detector. The electrical signal from 
the detector is analyzed and compared to the threshold value. A threshold voltage is 
defined to trigger the on/off performance of the solenoid 3 -way isolation valve 
installed behind the sensor, so as to allow the droplets of interest to be sorted into the 
corresponding channel. A hold-time and delay- time are set to determine how long the 
valve to be kept open and when the opening of the valve will be delayed after receiving 
a signal trigger from the GMR sensor. This feature is crucial when the spacing 
between droplets is smaller than the distance between the GMR sensor and the sorting 
junction. 

Multiparametric analysis. The detection events of ferrofluid droplets are represented 
by the detected peaks (Supplementary Fig. la). Each peak can be characterized by its 
amplitude, the full- width at half maximum (FWHM), positive/negative amplitude 
ratio (PNR). These characteristics can serve as fingerprints for each ferrofluid droplet. 
For instance, the amplitude is related to the size and concentration of droplets, the 
FWHM contains the information about the speed and size of droplets, while the 
orientation of magnetization of the ferrofluid droplets is relevant to the PNR. These 
distinct features can be used to discriminate each labeled magnetic object. We 
performed the following experiments at a fixed flow conditions. Therefore, 
modifications of the FWHM are mainly attributed to the different size of droplets. 
Based on the previous studies, the size and concentration of ferrofluid droplets can be 
controlled and deterministically correlated to the derived voltage change. 

We performed multiparametric analysis on the ferrofluid emulsion droplets. The 
amplitude and FWHM of the detection peak were derived by a Gaussian fit to the peak 
and taken for the established parametric analysis. Since the amplitude and FWHM are 
related to the droplet size and flowing speed, the distribution of the droplet size and 
amount of ferrofluids under fixed flow condition can be evaluated. FWHM is 
dependent not only on the size of droplets but also on their velocities. Therefore, the 
combined multiparametric analysis of the amplitude and FWHM can be used for the 
evaluation of the droplet size as well as the flow conditions. The multiparametric 
density plots for the ferrofluid droplets (with concentration of 15 mg/ml) with dif- 
ferent volumes are summarized in Figure 2a (in the manuscript). The plots collect 
over about 100 droplets. The broad distribution of the FWHM and amplitude (black 
cloud, Fig. 2a, in the manuscript) is consistent with the real time detection results 
(Supplementary Fig. la). Also the corresponding fitted histogram displayed in 
Figure 2c (in the manuscript) indicates a clear polydispersity of droplet size. The 
droplet volume increases with the increase of flow rate of ferrofluids. The size of 
ferrofluid droplets is smaller than 1 mm for the flow rate of ferrofluid up to 200 nl/s. 
Thus the sensor is still working within the dynamic range as marked by the blue 
region in Supplementary Figure lb. The amplitude of voltage change is sensitive to 
the change of droplet size. Thus, a narrow distribution of the data point cloud and the 
decrease standard deviation of the amplitude and FWHM in the fitted histograms 
(Fig 2c, in the manuscript) indicates a more uniform distribution of droplet size at 
larger volumes. The standard deviations are derived by a Gaussian fit to the histo- 
grams of amplitude and FWHM, respectively. The elliptical circles in Figure 2a,b in 
the manuscript are guide to the eyes. The elliptical circle is located at the center of the 
Gaussian fits to the histograms of the amplitude and FWHM, while the long and short 
axis of the elliptical circles take the half width values of the Gaussian fit to the 
amplitude and FWHM histograms. With droplet volume of 270 nl, the sensor is 
saturated. The variation of voltage change amplitude can reflect the concentration 
variation. However, in this case, the analysis shows that the standard deviation of 
amplitude is below 2 |j,V (see Supplementary Fig. 2i), which is smaller than the 
binning size of the histogram and close to the noise level (around 1.5 |aV). This 
indicates that concentration variation per droplet contributes to the change of 
amplitude to lesser extent in the experiments. The histogram of amplitude can clearly 
reflect the size distribution of droplets when the sensor is still working in the dynamic 
region. However, as mentioned, when the size of droplets is larger than the sensor 
width, the sensor response is saturated with even larger size of droplets. Nonetheless, 
from the FWHM, the information of droplet size can still be extracted, which can be 
seen from the fact that the magnitude of FWHM increases as the droplet volume 
increases from 120 nl to 270 nl .We next performed the multiparametric analysis of 
droplets with different concentrations which were produced under the same flow 
condition (Qon = 500 nl/s, Qferro = 300 nl/s). The clouds of detection characteristics 
for droplets with different concentrations naturally fall into the same FWHM region 
(Supplementary Fig. 2e-h) indicating the consistent, homogeneous distribution of 
the droplet size for each specified concentration. The standard deviation of FWHM 
and amplitude for each specified concentrations levels off (Supplementary Fig. 2j), 



which indicates a relative consistent and homogenous droplets produced under the 
same flow conditions for each specified concentration. However, due to the scaling of 
voltage change with the magnetic moment, the center of the data point cloud is clearly 
shifted to higher voltage amplitudes for droplets with larger concentration of mag- 
netic nanoparticles. 
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